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Abstract 

Simple closed-box (CB) models of chemical evolution are extended on 
two respects, namely (i) simple closed- (box+reservoir) (CBR) models 
allowing gas outflow from the box into the reservoir (Hartwick 1976) 
or gas inflow into the box from the reservoir (Caimmi 2007) with 
same composition as the preexisting gas and rate proportional to the 
star formation rate, and (ii) simple multistage closed- (box+reservoir) 
(MCBR) models allowing different stages of evolution characterized 
by different inflow or outflow rates. In any case, the stellar initial 
mass function is assumed to be universal, and mass conservation holds 
for the whole system (box+reservoir) while it is violated for each 
subsystem (box and reservoir). The theoretical differential oxygen 
abundance distribution (TDOD) predicted by the model, under the 
assumption of instantaneous recycling, is a continuous broken line, 
where outflow and moderate inflow rates are represented by negative 
slopes, steady inflow rates by null slopes, and strong inflow rates by 
positive slopes. Then an application is made to a special stellar sys- 
tem resembling the inner Galactic halo. To this aim (a) a fictitious 
sample is built up from two distinct samples (Ryan and Norris 1991; 
Schorck et al. 2009) and taken as representative of the inner Galactic 
halo, and (b) different [0/H]-[Fe/H] empirical relations are deduced 
from different samples (Rich and Boesgaard 2009; Fabbian et al. 2009; 
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Schmidt et al. 2009) related to different methods, and two of them are 
selected for determining the empirical differential oxygen abundance 
distribution (EDOD) with regard to the fictitious sample. In both 
cases the EDOD is represented, to an acceptable extent, as a con- 
tinuous broken line. The slopes and the intercepts of the regression 
lines are determined, and then used as input parameters to MCBR 
models with fiducial values assigned to the remaining input parame- 
ters. Output parameters are the gas mass fraction and the star mass 
fraction at the end of each evolutionary stage. If the inner halo and 
the metal-poor bulge (after the inflow stage) are represented by the 
box and the reservoir, respectively, then the inner halo fractional mass 
(normalized to the halo) is comparable with, or exceeds by a factor up 
to 4, the metal-poor bulge fractional mass (normalized to the bulge), 
for current estimates of the halo-to-bulge mass ratio of about 0.05- 
0.10. On the other hand, quantitative results cannot be considered 
for applications to the inner Galactic halo, unless selection effects and 
disk contamination are removed from halo samples, and discrepan- 
cies between different oxygen abundance determination methods are 
explained. 

keywords - galaxies: evolution - stars: formation; evolution. 

1 Introduction 

The empirical metallicity distributions of long-lived stars of different popu- 
lations, constrain models for the formation and the evolution of the Galaxy. 
Simple closed-box (CB) models make a useful tool in the description of galac- 
tic chemical evolution. The original formulation (Searle and Sargent 1972; 
Pagel and Patchett 1975) relies on instantaneous recycling and homogeneous 
mixing approximations i.e. gas from dying stars is instantaneously returned 
to and homogeneously mixed with the interstellar medium. In addition, (1) 
a null star mass fraction and a null metal abundance are taken as initial 
values, and (2) mass conservation is assumed to hold, which implies no gas 
outflow from the box or inflow into the box. 

Simple CB models may be extended in many respects, such as nonzero 
initial metal abundance (prompt initial enrichment) and/or star mass frac- 
tion (Truran and Cameron 1971), inhomogeneous mixing (Caimmi 2000 + 
erratum 2001a, hereafter quoted together as COO; Caimmi 2001b, hereafter 
quoted as C01), and gas outflow (Hartwick 1976, hereafter quoted as H76) 
or inflow (Caimmi 2007, hereafter quoted as C07) with same composition as 
the preexisting gas and rate proportional to the star formation rate. For- 
mulations implying more extreme changes, such as reject of instantaneous 
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recycling and unconditioned gas outflow or inflow, can no longer be consid- 
ered as "simple" models. 

The existence of a G-dwarf problem i.e. detection of too few metal defi- 
cient G dwarf (or, more generally, of a selected spectral type) with respect to 
that which could be expected from simple CB models of chemical evolution 
(e.g., Searle and Sargent 1972; Pagel and Patchett 1975; Haywood 2001) was 
first established in the solar neighbourhood (van den Bergh 1962; Schmidt 
1963). Although in a less extreme form, a G-dwarf problem appears to exist 
both in the halo (e.g., H76; Prantzos 2003) and in the bulge (e.g., Ferreras et 
al. 2003). In addition, a G-dwarf problem has been recognized in both bulge- 
dominated and disk-dominated galaxies (Henry and Worthey 1999), which 
is consistent with the idea that the G-dwarf problem is universal (Worthey 
et al. 1996). 

According to current ACDM cosmological scenarios, galaxies were largely 
built out of disrupted smaller subunits (dSph galaxies whose surviving cores 
could be massive globular clusters). A similar process is presently occurring 
on a larger scale: a central bulge (a cD galaxy) is accreting in galaxy clusters, 
at the expense of infalling smaller galaxies. In this view, a universal G- 
dwarf problem appears to be closely related to the initial assembling stage 
of galactic evolution. 

With regard to the Galaxy, the empirical differential oxygen abundance 
distribution (EDOD) shows that (i) at least two different stages of (chemical) 
evolution exist, and (ii) each stage can be described using an extended simple 
CB model, for both the inner halo (C07), the bulge (C07), the thick disk 
(Caimmi and Milanese 2009, hereafter quoted as CM09), and the thin disk 
(CM09). In particular, the earlier stage is marked by an initially increasing 
EDOD. 

The advantage of the EDOD, tp = log [AAy(NA0)@, on the oxygen abun- 
dance distribution, \& = AN/N, is that the theoretical differential oxygen 
abundance distribution (TDOD) predicted by extended simple CB models 
is a straight line on the (0<j)ip) plane, which provides a more sensitive test 
(Pagel 1989; Malinie et al. 1993; Rocha-Pinto and Maciel 1996; COO; C01; 
C07; Caimmi 2008; CM09). 

As shown in recent attempts (Carollo et al. 2007, 2010), the Galactic halo 
may be divided into two structural components, namely (1) an inner spheroid 
with pa 15 kpc, c ~ 10 kpc, which exhibits a modest prograde rotation 

and a metallicity peak at [Fe/H] pa —1.6, and (2) an outer sphere with radius, 

Caption of symbols: 4> = Zo/(Zo)q is the oxygen abundance normalized to the solar 
value, N is the population of the sample under consideration, and AiV is the number of 
sample objects within a bin centered in <fr and bounded at <p T A0/2. 
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R ps 20 kpc, which exhibits a clear retrograde net rotation and a metallicity 
peak at [Fe/H] ~ —2.2. In addition, the inner halo population dominates 
within its own volume, while the outer halo population includes a larger 
fraction of low-metallicity ([Fe/H] < —2.0) stars than does the inner halo 
population. 

For a component of orbital motion measured with respect to the Local 
Standard of Rest, V < 200kms _1 , the distribution of [Fe/H] appears similar 
to what in the past was considered "the halo" , with a single metallicity peak 
at [Fe/H] ps —1.6 (Carollo et al. 2007, Fig. 2), which implies the inner halo 
population is dominant in the related sample. A similar trend is exhibited by 
the metallicity distribution inferred from an earlier attempt (Ryan and Norris 
1991, hereafter quoted as RN91, Fig. 5c), and the related sample (hereafter 
quoted as the RN91 sample) can also be considered as representative of the 
inner halo. In both cases, data are biased towards low metallicities, [Fe/H] 
< -3.0. 

The inner halo low-metallicity tail has been determined using data from 
the Hamburg/ESO survey (Schorck et al. 2009, hereafter quoted as HV, 
Fig. 11), but the related sample (hereafter quoted as the H V sample) is bi- 
ased towards higher metallicities, [Fe/H] > —3.0, and extremely low metallic- 
ities, [Fe/H] < —4.2. In the comparison of the empirical with the theoretical 
metallicity distribution, selection effects must be taken into account. Ac- 
cordingly, each theoretical metallicity distribution under consideration has 
to be converted into its counterpart as it would be observed by applying the 
metal-poor star selection criteria used in the survey (HV). 

If, on the other hand, the HV sample is dominated by the inner halo 
population, a correction for the above mentioned selection effects is expected 
to yield a metallicity distribution with a single peak at [Fe/H] ps —1.6, similar 
to its counterpart related to earlier samples (RN91; Carollo et al. 2007, Fig. 2, 
V < 200 km s -1 ). In fact, the metallicity distributions related to the RN91 
and HV sample can be adjusted to match in the range, —3.4 < [Fe/H] < 
—2.6, where selection effects are negligible (HV, Figs. 11, 12). 

The Milky Way inner halo offers a unique opportunity to construct the 
star formation and mass assembly history of a galactic inner halo, hence 
providing a unique benchmark for theories of galaxy formation and evolution. 
To this aim, the sample used has to be complete at least for a wide range 
of metallicity. In the case under discussion, —3.6 < [Fe/H] < —2.8 for 
the HV sample and -3.4 < [Fe/H] < -1.0 for the RN91 sample. The 
remaining ranges are affected by various biases, mainly due to the occurrence 
of selection effects for the low-metallicity tail above a threshold (HV) and 
disk contamination for the high-metallicity tail (RN91, H V). 

The metallicity distribution of the inner halo can be deduced from a 
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fictitious sample (hereafter quoted as the fslO sample) within the range, 
—4.2 < [Fe/H] < +0.2, under a number of restrictive but reasonable as- 
sumptions, namely (i) the HV sample is representative of the inner halo 
within the range, —4.2 < [Fe/H] < —3.0; (ii) the RN91 sample is represen- 
tative of the inner halo within the range, —2.8 < [Fe/H] < +0.2; (iii) the 
H V and RN91 samples are equally representative (even if related to different 
populations) of the inner halo within the range, —3.0 < [Fe/H] < —2.8. 

In absence of direct [O/H] evaluations, the EDOD related to the fslO 
sample can be deduced by use of a [0/H]-[Fe/H] empirical relation which, 
on the other hand, is strongly dependent on both the selection of spectro- 
scopic oxygen lines and the choice of the model, as shown in recent attempts 
(Ramirez et al. 2007; Rich and Boesgaard 2009, hereafter quoted as RB09; 
Fabbian et al. 2009, hereafter quoted as Fal09; Schmidt et al. 2009, hereafter 
quoted as Sal09) where earlier attempts are quoted (for further insight refer 
to the proceedings edited by Barbuy et al. 2001). The discrepancy between 
[0/H]-[Fe/H] empirical relations, deduced using different methods and differ- 
ent models, still remains large. For this reason, the TDOD should be fitted 
to the EDOD related to two different [O/H]- [Fe/H] empirical relations, as 
done in previous papers (C01; C07; Caimmi 2008; CM09). 

Simple CB models of chemical evolution were first extended allowing for 
gas outflow (H76) and later for moderate gas inflow, with same metal abun- 
dance as in the preexisting gas, yielding negative TDOD slopes in the (O^ip) 
plane. On the other hand, the low-met allicity tail of the EDOD is fitted to a 
straight line with positive slope for both the halo (C07), the bulge (C07), the 
thick disk (CM09), and the thin disk (CM09), which is achieved by extending 
simple CB models to strong gas inflow with same metal abundance as the 
preexisiting gas. 

In this view, different stages of evolution are related to different domains 
in the normalized oxygen abundance, (0u)i < < (0u)/ ; where the EDOD 
is fitted to a regression line, ip = a\j(j) + by, with slope, a\j, and intercept, 
bv (U = I, II, is the stage considered, beginning and ending at i and / 
configurations, respectively). For reasons of continuity, the final values of 
parameters related to a selected stage must necessarily coincide with the 
corresponding initial values related to the next stage, with the exception of 
the outflow or inflow rate, which can be deduced from the EDOD. Accord- 
ingly, the knowledge of present values allows the calculation of past values, 
in the light of the model, getting insight on the formation and evolution of 
the system under consideration. 

In summary, the current paper is aimed to the following main points: (i) 
metallicity distribution related to a fictitious sample (fslO), supposed to be 
representative of the inner halo (subsection I2.ip ; (ii) choice of two different 
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[0/H]-[Fe/H] empirical relations deduced from recent samples (RB09, Fal 
09, Sal09) (subsection 12. 2p ; (iii) EDOD determination from the fslO sample 
by use of the above mentioned [0/H]-[Fe/H] empirical relations (subsection 
12. 3ft : (iv) formulation of extended simple CB models, namely (a) simple CBR 
models, allowing for gas outflow or inflow with same metal abundance as the 
preexisting gas, and (b) simple MCBR models, allowing for different stages 
of evolution (subsections 13. 1) 13.21 and 13.4)) : (v) determination of the best 
fitting TDOD to the EDOD inferred from the fslO sample (subsection 13. 3ft : 
(vi) application to a special stellar system resembling the inner Galactic halo 
(subsection 13.5)1 . In addition, the discussion and the conclusion make the 
subject of subsection 13.61 and section H], respectively. 

2 Data and inferred quantities 

The following ingredients are necessary for the determination of the EDOD 
(in particular, related to the inner Galactic halo), from which input parame- 
ters of a selected model are deduced: (i) a representative sample; (ii) a related 
[Fe/H] distribution; (iii) a [0/H]-[Fe/H] empirical relation. Each point shall 
be separately dealt with in the forthcoming subsections. 

2.1 Building up a fictitious sample (fslO) 

In building up a fictitious sample (the fslO sample) of inner halo stars, two 
different samples shall be considered, namely the RN91 sample (N = 372) 
of kinematically selected halo subdwarfs and the HV sample (N = 3439) 
of metal-poor stars selected from the Hamburg/ESO objective prism survey. 
More specifically, the HV sample has been determined from a subsample 
(N = 1638) with available spectroscopic follow-up observations, by means of 
scaling to the full parent sample (N = 3439). For further details refer to the 
parent paper (HV). The RN91 sample can be considered as complete within 
the range, —3.4 < [Fe/H] < —1.0. At different abundance the sample is 
incomplete and suffers contamination from disk stars at the high-metallicity 
tail, [Fe/H] > —1.0. The HV sample can be considered as complete within 
the range, —4.2 < [Fe/H] < —2.8. At different abundance the sample is 
incomplete and suffers both selection effects ([Fe/H] > —2.8) and contami- 
nation from disk stars ([Fe/H] > —2.0). 

Contamination from outer halo stars also exists in both samples, but the 
effect is expected to be negligible. In fact, the inner halo dominates within 
a spheroid centered on the Galactic centre, with equatorial plane coinciding 
with the Galactic plane, and semiaxes (a, c) = (15, 10) kpc (Carollo et al. 
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2007). On the other hand, both the RN91 and H V sample are made mainly of 
stars placed within the above mentioned spheroid. Then the contamination 
from outer halo stars is expected to be appreciable only at low ([Fe/H] < 
—2.0) metal abundances (Carollo et al. 2007). 

The metallicity distribution from the RN91 sample peaks at [Fe/H] w 
— 1.6, in agreement with recent results for the inner halo population (Car- 
ollo et al. 2007, 2010). On the contrary, the metallicity distribution related 
to the HV sample appears to be bimodal with an absolute maximum at 
[Fe/H] w —2.2 and a relative maximum at [Fe/H] « —1.0. The bimodality 
is probably an artefact due to selection effects and/or disk contamination 
for the following reasons^): (i) the HV sample is mainly made of inner halo 
stars (HV); (ii) the expected peak at [Fe/H] w —1.6 (Carollo et al. 2007, 
2010) is lacking; (iii) the metallicity distribution from the RN91 sample can 
be scaled to match its counterpart from the HV sample within the range, 
—3.4 < [Fe/H] < —2.6 (HV); (iv) a bimodal metallicity distribution cannot 
be fitted to its theoretical counterpart deduced from simple models of chem- 
ical evolution, even if the latter is "observed" by applying the metal-poor 
selection criteria used in the Hamburg/ESO survey (HV). 

The H V sample is mainly made of giant stars. A subsequent study on 
the stellar content of the Hamburg/ESO survey has been focused on a sam- 
ple (N = 617) of main sequence turnoff stars (Li et al. 2010). Both samples 
exhibit a quite analogous metal abundance distribution where the halo popu- 
lation dominates ([Fe/H] < —2.0), while the contrary holds for higher values 
([Fe/H] > —2.0). In the latter case, the H V sample is also affected from the 
survey- volume effect. For further details refer to the parent paper (Li et al. 
2010). 

Basing on the above considerations, the following working hypotheses are 
made. 

(1) The H V sample (N = 3439) is representative of the inner halo within the 

range, —4.2 < [ Fe/H] < —2.8, where selection effects are minimized 
and contamination from disk stars is negligible (HV). 

(2) The RN91 sample (N = 372) is representative of the inner halo within 

the range, —3.0 < [ Fe/H]< +0.2, with a caveat due to contamination 
from disk stars for [Fe/H] > -1.0 (RN91). 

(3) The RN91 and HV sample are equally representative (even if belonging 

2 The problem was raised via e-mail (December 2009) to two coauthors of the parent 
paper (HV) and to two coauthors of the paper where inner and outer halo populations 
were first detected (Carollo et al. 2007), but no reply was received up to-day. It cannot 
be excluded that related messages have been lost into the spam bin. 
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to different populations) of the inner halo within the range, —3.0 < 
[Fe/H] < —2.8, where the number of stars is (AiV) RN91 = 8 and 
(AiV) H v = 160, respectively. 

Accordingly, the fslO sample can be built up from the RN91 and H V samples, 
by normalizing to the same number of stars within the metallicity range where 
the above mentioned samples are supposed to be equally representative of 
the inner halo. More specifically, the normalization shall be performed with 
respect to the H V sample, but the Poissonian error related to each bin shall 
remain unchanged with respect to the parent sample, with the minimum 
among the two retained for the bin where the parent samples are equally 
representative of the inner halo. Then the number of stars belonging to the 
fslO sample, with regard to a selected metallicity bin, remains unchanged or 
is multiplied by a factor 20 according if the parent sample is HV ([Fe/H] < 
-2.8) or RN91 ([Fe/H] > -2.8), yielding a fictitious population, N = 7452. 

2.2 The [0/H]-[Fe/H] empirical relation 

In dealing with simple models of chemical evolution, involving the assump- 
tion of instantaneous recycling, the predicted metal abundance has to be 
compared to the observed oxygen (or any other primary element) abundance 
(e.g., Pagel 1989; COO; C01). Unfortunately, oxygen is more difficult than 
iron to detect, and an empirical formula may be needed to express the former 
as a function of the latter. The population of available samples where oxygen 
abundances are directly determined, does not exceed a few hundreds at most 
(e.g., Ramirez et al. 2007; Melendez et al. 2008; RB09; Fal09; Sal09). With 
regard to the halo, only the RB09 and Fal09 samples can be considered as 
representative, while the Sal09 sample shall be included for comparison, and 
the remaining above quoted two shall be excluded. 

The RB09 sample (N = 49) is made of a homogeneous subsample (N = 
24) of metal-poor (—3.5 < [Fe/H] < —2.2) stars, and a non homogeneous 
subsample (N = 25) of higher-metallicity (—3.1 < [Fe/H] < —0.5) stars. 
In both cases, the stellar population remains unspecified and oxygen abun- 
dance has been determined using standard local thermodynamical equilib- 
rium (LTE) one-dimensional hydrostatic model atmospheres. The calculated 
oxygen abundance is independent of the solar value. Typical uncertainties 
are A [Fe/H] = =f0.15 and A[0/H] = =(=0.15. Standard deviations are also 
provided for each star. For further details refer to the parent paper (RB09). 

The Fal09 sample (N = 43) is made of halo stars (-3.3 < [Fe/H] < -1.0) 
where oxygen abundance has been determined using three different meth- 
ods involving (a) LTE one-dimensional hydrostatic model atmospheres; (b) 



8 



Table 1: Values of effective temperature, T eff , decimal logarithm of surface 
gravity, logg, iron abundance, [Fe/H], and oxygen abundance, [0/H], for 
N = 11 stars in common between the RB09 and Fal09 sample. With regard 
to the latter, oxygen abundance has been determined using three different 
methods, LTE, SHO, SHI, respectively, and oxygen abundances related to 
LP831-070 are upper limits. For further details refer to the text. 



star 


T c s 


RN09 
logs [Fe/H] 


[O/H] 
LTE 




logfi- 


Fal09 
[Fe/H] [O/H] 
LTE 


[O/H] 
SHO 


[O/H] 
SHI 


LP651-04 


6030 


4.26 


-2.89 


-2.04 


6371 


4.20 


-2.63 


-1.62 


-2.21 


-1.93 


BD-13°3442 


6090 


4.11 


-2.91 


-2.15 


6366 


3.99 


-2.69 


-1.77 


-2.39 


-2.11 


G 11-44 


5820 


3.58 


-2.29 


-1.63 


6178 


4.35 


-2.03 


-1.29 


-1.47 


-1.40 


G64-12 


6074 


3.72 


-3.45 


-2.24 


6435 


4.26 


-3.24 


-2.21 


-3.10 


-2.71 


G64-37 


6122 


3.87 


-3.28 


-2.32 


6432 


4.24 


-3.08 


-2.24 


-3.12 


-2.70 


LP635-14 


5932 


3.57 


-2.71 


-2.00 


6367 


4.11 


-2.39 


-1.60 


-2.03 


-1.85 


LP815-43 


6405 


4.37 


-2.76 


-1.86 


6483 


4.21 


-2.71 


-1.95 


-2.70 


-2.36 


HD084937 


6206 


3.89 


-2.20 


-1.49 


6357 


4.07 


-2.11 


-1.39 


-1.64 


-1.56 


HD140283 


5692 


3.47 


-2.56 


-1.72 


5849 


3.72 


-2.38 


-1.67 


-1.91 


-1.81 


HD194598 


5875 


4.20 


-1.25 


-1.00 


6020 


4.30 


-1.15 


-0.51 


-0.75 


-0.68 


LP831-070 


5985 


4.75 


-3.06 


-2.45 


6232 


4.36 


-2.93 


-2.13 


-2.95 


-2.54 



three-dimensional hydrostatic model atmospheres in absence of LTE with no 
account taken of the inelastic collisions via neutral H atoms (Su = 0), here- 
after quoted as SHO; (c) three-dimensional hydrostatic model atmospheres in 
absence of LTE with due account taken of the inelastic collisions via neutral 
H atoms (Su = 1), hereafter quoted as SHI. For a single object (LP831-070) 
only an upper limit to oxygen abundance has been determined. Typical un- 
certanties are A [Fe/H] = =1=0.15 and A [O/H] = =f=0. 15. Standard deviations 
are not provided for each star. For further details refer to the parent paper 
(Fal09). 

The RB09 and Fal09 samples have N — 11 (necessarily halo) stars in com- 
mon, where the values assumed for effective temperature and surface gravity 
have been determined using different methods, yielding different values for 
each star, as shown in Table[TJ A comparison between [0/H]-[Fe/H] empirical 
relations deduced from the data listed in Table [TJ, in connection with RB09 
and Fal09 (case LTE, left; SHI, right) samples, is presented in Fig.[TJ, upper 
panels, while the correspondance between [Fe/H] (left) and [O/H] (right), 
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deduced from the above mentioned data, is shown in lower panels. With 
regard to Fal09 sample, case SHO has not been considered, as it seems to 
overcorrect LTE abundances and yield values of [O/Fe] which are probably 
too low. For further details refer to the parent paper (Fal09). 

It can be seen that the [0/H]-[Fe/H] empirical relation is slightly affected 
by the parent sample in the LTE case, with the exception of the most metal- 
rich sample object (Fig.[TJ top left panel). A low discrepancy appears in the 
SHI case (Fig.HJ top right panel), where results in absence of LTE approxima- 
tion are available for the Fal09 sample only. The occurrence of a systematic 
difference in [Fe/H] determinations is probably due to different methods re- 
lated to different samples, which makes all points lie above the straight line 
with unit slope (Fig.[U bottom left panel). The same holds, though to a 
lesser extent, for [O/H] determinations related to LTE case (Fig.[[J, bottom 
right panel, squares), while points related to the SHI case (Fig.[TJ bottom 
right panel, crosses) lie both above and below the straight line with unit 
slope. 

The normalized [Fe/H] distributions related to the RB09 and Fal09 sam- 
ple are plotted in Fig. [2] (upper panels), where a similar trend is shown and, 
in particular, both exhibit a peak near [Fe/H] = —2.2, similar to their coun- 
terpart deduced from the HV sample (bottom right panel) but in contrast 
with their counterpart deduced from the RN91 sample (bottom left panel). 
Then the RB09 and Fal09 samples seem to be more representative of the 
outer halo than the inner halo, where the most populated metallicity bin 
is placed at [Fe/H] —2.2 and [Fe/H] rj —1.6, respectively (Carollo et al. 
2007, 2010). Accordingly, the determination of the EDOD from the [Fe/H] 
distribution related to RN91 and HV samples has to be made under the as- 
sumption that the [0/H]-[Fe/H] empirical relations, deduced from the RB09 
and Fal09 samples, hold to a similar extent for both the inner and the outer 
halo. 

The following regression models are used for fitting to [O/H]- [Fe/H] em- 
pirical relations. 

G heteroscedastic functional models where the errors in X and in Y are 
uncorrelated (York 1966). 

Y homoscedastic structural models where the errors in X are negligi- 
ble (ideally null) with respect to the errors in Y (Isobe et al. 1990, 
OLS(Y|X) thereiii). 

3 The captions, OLS(Y|X) and OLS(X|Y), have to be interchanged one with the other 
in Table 1 of the quoted paper. 
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X homoscedastic structural models where the errors in Y are negligi- 
ble (ideally null) with respect to the errors in X (Isobe et al. 1990, 
OLS(X|Y) therein?). 

B homoschedastic structural models where the regression line bisects the 
angled formed by the regression lines related to models Y and X above 
(Isobe et al. 1990, OLS bisector therein). 

O homoscedastic structural or functional models where the regression line 
minimizes the sum of the perpendicular distances between the data 
points and the line (Isobe et al. 1990; Feigelson and Babu 1992; or- 
thogonal regression therein). 

R homoscedastic structural models where the regression line has slope 
equal to the geometric mean of the slopes of the regression lines re- 
lated to models Y and X above (Isobe et al. 1990, reduced major-axis 
therein) . 

The [0/H]-[Fe/H] empirical relations deduced from RB09, Fal09 (cases 
LTE, SHO, SHI), and Sal09 samples, shown in Fig.[3j are interpolated using 
the regression models listed above. For heteroscedastic data and homoscedas- 
tic models, the typical uncertainties are assigned to all the data points. The 
results are listed in Table [2j An inspection of Fig. [3] and Table 2 shows that 
systematic errors related to different methods and/or different spectral lines 
in oxygen abundance determination, are dominant on both the dispersion 
due to measurement errors and intrinsic scatter. At this stage, a precise fit 
to the data extracted from a selected sample would be meaningless, and only 
acceptable fits related to different situations shall be considered. Accordingly, 
the following [0/H]-[Fe/H] empirical relations: 

[O/H] = 0.72[Fe/H] ; (1) 

with regard to the RB09 sample, Fig. [3] (full lines on data points), and 

[O/H] = [Fc/H] + 0.70 ; (2) 

with regard to the Fal09 sample, case SHI, Fig.[3] (dashed lines on data 
points), shall be used for determining the EDOD from the [Fe/H] distribu- 
tion, related to a selected sample. The regression lines related to five different 
methods listed in Table El are shown for each sample on the top right panel 
ofFig.EJ 

4 The captions, OLS(X|Y) and OLS(Y|X), have to be interchanged one with the other 
in Table 1 of the quoted paper. 
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Table 2: Regression line slope and intercept estimators, a and b, and related 
dispersion estimators, <3"a, and <7g, for the models (m) mentioned in the text, 
applied to the [0/H]-[Fe/H] empirical relation deduced from the following 
samples: RB09 (top left); Fal09, cases LTE (top right), SHO (middle right), 
SHI (bottom right); and Sal09 (middle left). 



m a 6"a b a~ b a b a b 
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0665 
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,0310 
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0706 
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0037 
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+0. 


.1431 
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0. 
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1. 
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0. 


.0441 


0. 
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0. 
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0. 
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0. 
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0. 
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0. 
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2.3 The empirical differential oxygen abundance dis- 
tribution (EDOD) 

With regard to a selected [0/H]-[Fe/H] empirical relation: 

[0/H] = a[Fe/H] + b ; (3) 

and a specified [Fe/H] distribution, let [Fe/H], [Fe/H] T , be the coordinates of 
a selected bin centre and bin left (— ) and right (+) extremum, respectively, 
and [0/H], [0/H] T their counterparts deduced via Eq. <^5§ for the related 
[0/H] distribution. Accordingly, the bin semiamplitude reads: 

AT[0/H ] = 1°/ H 1 + - 1°/"!' =a [Fe/H]^-[Fe/Hl- =aAT|Fg/H] . ^ 

where [Fe/H] and [0/H] are logarithmic number abundances normalized to 
the solar value (e.g., C07). 

The oxygen mass abundance normalized to the solar value, <ft, to a good 
extent may be expressed as (e.g., C07): 

Z ° exp 10 [O/H] ; (5) 



r (Z Q )e 

and a selected bin centre and bin semiamplitude, related to the [O/H] dis- 
tribution, read (e.g., C07): 

= l{exp 10 [0/H] + + exp 10 [0/H]-} ; (6) 
A^</> = i{ex Pl0 [0/H] + - ex Pl0 [0/H]-} ; (7) 

where the bin width is variable for a constant bin width related to the [0 /H] 
distribution. The [0/H]- [Fe/H] and 0-[Fe/H] empirical relations, expressed 
by Eqs. (pQ) and fl2]), by use of Eqs. flB]) and ([7j) in the latter case, are rep- 
resented in Table [31 where Bp = [Fe/H] and -Bo = [0/H] to save space. 

Following recent attempts (Rocha-Pinto and Maciel 1996; COO; C01; 
C07), the EDOD in a selected class of objects is defined as: 

AN 

i>{4>) = lo ; 



NA<p ' v ' 

where AiV is the number of objects within a normalized oxygen abundance 
bin, A(f), centered in <ft, and iV is the number of sample objects. The increment 
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Table 3: The normalized oxygen abundance, = Zo/{Zo)&, for two different 
[0/H]-[Fe/H] empirical relations deduced from interpolation to two different 
data sets, namely RB09 in presence of the local thermodynamic equilibrium 
approximation (left side), and Fal09 in absence of the local thermodynamic 
equilibrium approximation, case SHI (right side). Upper and lower values for 
each bin are denoted by the apex + and — , respectively. The corresponding 
bin mean value and semiamplitude, with regard to normalized oxygen abun- 
dance, are labelled as <fi and A T 0, respectively. Labels B^: and B^ stand for 
[Fe/H] T and [0/H] T , respectively. 
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ratio, AN/A(p, used in earlier attempts (Pagel 1989; Malinie et al. 1993) 
is replaced by its normalized counterpart, AN/(NA<f)), to allow comparison 
between different samples. The uncertainty on AiV is evaluated from Poisson 
errors (e.g., RN91), as <taat = (AN) 1 / 2 , and the related uncertainty in the 
EDOD is: 



A^ = \ip — = log 



(AiV) 1/2 

1 =f — 

^ AiV 



(9a) 



-p , AiV =f (AiV) 1 / 2 

where — >■ — oo in the limit AiV — > 1. For further details refer to the 
parent papers (C01; C07). 

The EDOD related to RN91, HV, and fslO samples are listed in Table 
H]for the [0/H]-[Fe/H] empirical relations expressed by Eqs. ([T]) and ([2]), left 
and right side, respectively. The EDOD related to the fslO sample, taken to 
be representative of the inner halo, is plotted in Fig.[4]for the [0/H]-[Fe/H] 
empirical relations expressed by Eqs. (DP and ([2]), left and right panel, re- 
spectively. Upper panels represent the whole distribution, while lower panels 
zoom the low-metallicity range. The vertical bars on the horisontal axis 
mark [Fe/H] = —2.2, —1.6, where the global [Fe/H] distribution of the outer 
and inner halo, respectively, peak according to recent results (Carollo et al. 
2007, 2010), and [Fe/H] = —0.8, where a transition from halo to bulge/disk 
globular clusters occurs (Mackey and van den Bergh 2005). 

A main feature of the EDODs plotted in Fig.HJ is the presence of five 
regions characterized by a nearly linear trend, which shall be named O, A, 
F, C, E, respectively, and defined in the following metallicity ranges, each 
containing nx bins, X=0, A, F, C, E. 

O [Fe/H] < -4.2; [Fe/H] > -0.2; n -> +oo; the distribution coincides 
with the horisontal axis after removing a single high-metallicity star 
from the RN91 sample, considered as due to disk contamination or, in 
any case, an outlier. On the other hand, the last appears in related 
tables and figures for comparison. 

A —4.2 < [Fe/H]< —2.7; n\ = 6; the distribution is steep with positive 
slope. 

F —2.7 < [Fe/H] < —1.7; n F = 6; the distribution is mild with negative 
slope. 

C —1.7 < [Fe/H]< —0.9; uq = 5; the distribution is steep with negative 
slope. 
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Table 4: The empirical, differential oxygen abundance distribution (EDOD) 
in the inner halo, deduced from the fslO sample (N = 7452) using two differ- 
ent [0/H]-[Fe/H] empirical relations, determined from interpolation to two 
different data sets, RB09 in presence of the local thermodynamic equilibrium 
approximation (left side), and Fal09 in absence of the local thermodynamic 
equilibrium approximation, case SHI (right side). The fictitious fslO sam- 
ple results from the combination of the HV sample (N = 3439) for lower 
metallicities, -4.2 <[Fe/H]< -3.0, and the RN91 sample (N = 372) for 
higher metallicities, —2.8 <[Fe/H]< +0.2, under the assumption that the 
two samples are equally representative of the inner halo within the metallic- 
ity bin, —3.0 < [Fe/H]< —2.8. The error on the generic bin height has been 
estimated from the Poissonian error of its counterpart related to the parent 
sample. See text for further details. 
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E —0.9 < [Fe/H] < —0.2; n E = 4; the distribution is less steep with 
negative slope. 

In absence of LTE approximation with regard to the [0/H]-[Fe/H] empirical 
relation deduced from the Fal09 sample (case SHI), regions F and C may 
merge into a single region, FC, within the range, —2.7 < [Fe/H] < —0.9, 
including npc — 10 bins. The vertical dashed lines in Fig. H] mark [Fe/H] = 
—2.7, —1.7, —0.8, from the left to the right. It can be seen that the global 
[Fe/H] distribution peaks early within the F and C region for the outer and 
inner halo, respectively. In addition, data points on the boundary between 
adjacent regions follow the linear trend exhibited by every of them. 

The regression line related to the EDOD within each populated region, 
has been determined using the B model listed in Table El under the assump- 
tion that the intrinsic scatter is dominant (Isobe et al. 1990). A single 
high-metallicity star from the RN91 sample, considered as due to disk con- 
tamination or, in any case, an outlier, has not been included in the fitting 
procedure. The regression line slope and intercept estimators and related 
dispersion estimators are listed in Table [5] for each region of the EDODs 
plotted in Fig.Hl The results are consistent (within =pr) with their counter- 
parts determined using the other models listed in Table [2J The selection of 
a special method is of little relevance, due to the paucity of data within each 
region. 

The regression lines are represented in Fig. [5] for each region (from the left 
to the right): A (dotted, positive slope), F(dotted, negative slope), C (full), 
E (dashed), with regard to the EDOD plotted in Fig.Hl In absence of LTE 
approximation, case SHI (right panels), a more inclined dashed line fits to 
the FC region, in alternative to F and C regions separately. To ensure con- 
tinuity, the above mentioned regions must be redefined by the intersections 
of regression lines, which make the transition from a selected region to the 
adjacent one. The results are listed in Table [6] for the EDOD and related 
regression lines plotted in Fig.0 A vertical line instead of a regression line 
has been considered for the O region, where no data exist. 

In conclusion, the EDOD related to the inner halo may be approximated, 
to a satisfactory extent, as the sum of four (A, F, C, E) or three (A, FC, 
E) regions, within which the trend is linear. An interpretations in terms of 
simple models of chemical evolution is highly attractive, as the corresponding 
TDOD shows, in fact, a linear trend (Pagel 1989; Rocha-Pinto and Maciel 
1996; COO; C01; C07). 
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Table 5: Regression line slope and intercept estimators, Ob and 6b, and re- 
lated dispersion estimators, <7a B , and <7g , for the B model applied to the 
oxygen abundance distribution (EDOD) plotted in Fig.|4]with regard to dif- 
ferent [0/H]-[Fe/H] empirical relations, deduced from the RB09 sample in 
presence of LTE approximation via Eq. flTJ (top panel), and from the Fal09 
sample in absence of LTE approximation, case SHI, via Eq. ([2]) (bottom 
panel). The method has been applied to each region (X) separately. Data 
points on the boundary between adjacent regions are used for determining 
regression lines within both of them. 
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FC 


-2.8643 E+0 


1.3796 E-l 


+8.2804 E-l 


2.7988 E- 


-2 


E 


-5.3661 E-l 


3.5681 E-2 


-6.1514 E-l 


4.3965 E- 


-2 



18 



Table 6: Transition points between adjacent regions, as determined from the 
intersection of related regression lines, for the oxygen abundance distribution 
(EDOD) plotted in Fig. |4] with regard to different [0/H]-[Fe/H] empirical re- 
lations, deduced from the RB09 sample in presence of LTE approximation 
via Eq. (pQ) (left panel) and from the Fal09 sample in absence of LTE approx- 
imation, case SHI, via Eq. (J2J) (right panel). In the latter case, the FC region 
has also been considered in alternative to F and C regions separately. For 
further details refer to the text. 





RB09 


Fal09 (SHI) 


transition 






O-A 


9.4624 E- 


-4 


+2.0383 E- 


3 


3.1623 E-4 


+3.0083 E- 


1 


A-F 


9.7001 E- 


-3 


+9.9844 E- 


1 


7.4271 E-3 


+8.6567 E- 


1 


F-C 


6.1687 E- 


-2 


+8.8953 E- 


1 


1.2740 E-l 


+9.9952 E- 


1 


C-E 


2.2854 E- 


-1 


-3.3795 E- 


-1 


6.4083 E-l 


-9.5902 E- 


1 


E-0 


7.1779 E- 


-1 


-1.4421 E- 


-0 


3.1623 E-0 


-2.3121 E- 





A-FC 










6.7114 E-3 


+8.0882 E- 


1 


FC-E 










6.2002 E-l 


-9.4785 E- 


1 



3 The model 

In their original formulation, simple models of chemical evolution are closed- 
box (CB) i.e. mass conservation (gas + stars) holds (e.g., Searle and Sargent 
1972; Pagel and Patchett 1975). In later formulation, a "reservoir" is added 
to the "box" , where mass conservation no longer occurs within the box, but 
still holds within the system box + reservoir. Accordingly, related models can 
be conceived as closed- (box+reservoir) (CBR) models. The gas within the 
box is "active" in the sense that allows star formation according to a specified 
birth-rate stellar function. The gas within the reservoir is "inhibited" in the 
sense that no star formation takes place. Gas may outflow from the box 
into the reservoir (H76) or inflow into the box from the reservoir (C07). 
The related TDOD is a straight line where the slope (to the knowledge of 
the author) is negative in all cases studied in literature (e.g., Pagel 1989; 
Rocha-Pinto and Maciel 1996; COO; C01; C07). 

For this reason, the formulation of CBR models shall be extended to a 
TDOD with positive slope, under the standard assumptions of CB models: 
(i) instantaneous recycling within the box, where stars are divided into two 
categories, namely (a) short-lived, which instantaneously evolve, and (b) 
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long-lived, all of which are still evolving; (ii) instantaneous mixing within 
the box, where the gas returned from short-lived stars is instantaneously 
mixed with the interstellar medium yielding uniform composition; and the 
standard assumptions of CBR models: (iii) mass conservation within the 
system (box+reservoir); (iv) gas outflow from the box into the reservoir or 
inflow into the box from the reservoir, at a rate proportional to the star 
formation rate; (v) inhibition of star formation within the reservoir; (vi) gas 
outflow or inflow with same composition with respect to preexisting gas. 

In this picture, the oxygen (or any other primary element) mass fraction 
can be determined, extending the procedure followed for CB models (e.g., 
Pagel and Patchett 1975; Wang and Silk 1993; COO) to CBR models. 



3.1 Basic theory 

The change in oxygen (or any other primary element) gas mass, M g o, is 
owing to four contributions, as: 

dM g0 _ / dMgcA / dM g0 \ / dM gC / 




dt \ dt I . I d* 7 V dt J a 

\ / si \ / gr \ / sdu 

related to subtraction via star formation (sf), subtraction via outflow from 
the box into the reservoir or addition via inflow into the box from the reservoir 
(gr), addition via unsynthesised gas from short-lived stars (sdu), and addition 
via newly synthesised gas from short-lived stars (sds), respectively. 
For simple CBR models, the following relations hold: 

dM g dM s dM s . d.Mg 

-a— — -na— — = -a(l + K)-^- ; (11) 

= - Zo a(l + K)^ ; (12) 

sdu 




(13) 

where M g is the mass in active gas, M s is the global mass in gas which 
has been turned into stars, M g o is the oxygen mass in active gas, a is the 
fraction in long-lived stars and stellar remnants within a star generation (lock 
parameter), k is the ratio of flowing (outflow from the box into the reservoir 
or inflow into the box from the reservoir) gas rate to locking (in the form 
of long-lived stars and stellar remnants) gas rate (flow parameter), Z Q and 
Z are the oxygen and metal mass abundance, respectively, within the active 
gas, and p is the ratio of the oxygen mass newly synthesised and returned 
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to the interstellar medium, for a metal-free initial composition^], to the mass 
locked up in long-lived stars and stellar remnants (yield parameter). The 
lock parameter, a, the flow parameter, k, and the yield parameter, p, for 
sake of brevity, in the following shall be quoted as the lock, the flow, and the 
yield, respectively. 

The definition of oxygen mass abundance, Zq, fractional active gas mass, 
H, and fractional star mass, s, read: 



M 



go 



M g 
Mo 



aS 



"MrT 



Ms 
M 



(14) 



where M s = aMs is the mass in long-lived stars and stellar remnants, and 
Mo is the total mass within the box at the starting configuration. 
The substitution of Eqs. (fTTT)- (fI5|) into (HE) yields: 



dM, 



go 



dt 



Zn- 



1 



1 + K 



P 



dt 



and the fractional oxygen gas mass, via Eq. (I14p reads: 

M g0 M g0 M g 



M Mg M 



Zqh ; 



(15) 



(16) 



where, in addition: 



dQgo/f) 
dt 



v d// dZ Q 



(17) 



on the other hand, in terms of fractional masses, Eq. (j!5p may be cast under 
the equivalent form: 



dQgo/x) 
dt 



-P 



d/x 
~dl 



and the combination of Eqs. (TT7|) and (|T8|) yields: 



o 



1 - Z 



,„d/i 



-P 



P 



P 



1 



K 



(19) 
(20) 



where p" is the effective yield parameter, hereafter quoted as the effective 
yield (H76). 



3 This detail is usually omitted in literature. 
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In the special case of a linear Zq-Z relation: 

A < 1 



1-Zq = 1 
l-Z A 



the combination of Eqs. (fT9|) and (12T]) produces: 

dZ Q __ Ap ,^ . 



1 - Zq fJL 

which can be integrated. After some algebra, the result is: 

Ap" 



(21) 



(22) 



Z - (Z Q )i = [1 - (Z Q )i 



(23) 



where the index, i, denotes the starting configuration at the cosmic time, t^. 
Reversing the role of the variables, Eq. (123 j) reads: 



E 

Pi 



1-Zq 
1 - (Z Q )i 



l/(Ap") 



(24) 



which is monotonic in Zq within the domain, < Zq < 1. 
Using the MacLaurin series development: 



E 



Ap" 



exp Ap" In 



1 + Af\n^ + 

Hi 



(25) 



under the further assumption that the terms of higher order with respect to 
the first can be neglected, Eq. (123]) reduces to: 





)i = 


-[1 


-( 


Zo) 


Ap"\nE 


< 


p" In 


E 


< 


Hi 









Hi 



p"\n^ 

Hi 



(26) 
(27) 



where the index, /, denotes the ending configuration at the cosmic time, tf, 
and Hi < H < Hf OT Hi > H > Hf owing to Eq. (|24|) . The initial oxygen 
and metal abundance may safely be neglected with respect to unity, (Zq)i < 
Zi 1, which implies A—>1 via Eq. (|2~TI) . Accordingly, Eq. (12^1) reduces to: 



Z - (Z )i = -p"\n— ; 

Hi 



(2? 



which is the classical formulation (H76). 
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It is worth emphasyzing that the inequality, expressed by Eq. (J27J), does 
not affect the instantaneous recycling approximation, but only the general 
formulation expressed by Eq. (1231) provided the Zq-Z relation is linear, ac- 
cording to Eq. ( 121]) . In general, it is assumed the instantaneou recycling 
approximation holds for sufficiently high fractional active gas mass fraction, 
/i ~ 0.1 (e.g., Prantzos 2007, Fig. 12 therein), but the threshold could be 
lowered using Eq. (f2"3~j) instead of Eq. (1251) . In fact, neglecting the terms of 
higher order in the MacLaurin series, expressed by Eq. (1251) . makes oxygen 
abundance, Zo, increased with respect to the general case, Eq. (1231) . which, 
in turn, is an upper limit with respect to the exact formulation, Eq. ffT9]) . 
More specifically, /i — > implies Z Q — >• 1 according to Eq. (123]) instead of 
Zq — > Zq({i = 0) < 1 due to the presence of other metals. Then the in- 
stantaneous recycling approximation expressed by Eq. ( T2"3~j) . or better by the 
integral of Eq. (|T9|) . appears to hold well even below the threshold, fi ~ 0.1. 

In terms of fractional masses, Eq. (fTT|) via Eq. ( TH]) may be cast under the 
equivalent form: 

. dS . ds du , , 

a( l + »)_ = (! + »)_ = --£ ; (29) 

which can be read in the following way: any mass change in active gas is 
counterbalanced by a change in long-lived stars and stellar remnants plus a 
change in gas outflow from the box into the reservoir {k > 0) or a change 
in gas inflow into the box from the reservoir (k < 0). More specifically, a 
number of different flow regimes may be distinguished as listed below. 

• Outflow regime (k > 0), where star formation is inhibited (H76). For 
an exhaustive description refer to earlier attempts (COO; C01). 

• Stagnation regime (n = 0), where star formation is neither inhibited 
nor enhanced. Accordingly, CBR models reduce to CB models. For 
an exhaustive description refer to earlier attempts (Searle and Sargent 
1972; Pagel and Patchett 1975). 

• Moderate inflow regime (—1 < k < 0), where star formation is en- 
hanced and active gas mass fraction monotonically decreases in time. 
For an exhaustive description refer to an earlier attempt (C07). 

• Steady inflow regime (k = —1), where star formation is enhanced and 
active gas mass fraction remains unchanged. 

• Strong inflow regime (k < —1), where star formation is enhanced and 
active gas mass fraction monotonically increases in time. 
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The effective yield, p", defined by Eq. (IZOl . cannot exceed the real yield, 
p, both in outflow and in strong inflow regime, while the contrary holds 
in moderate inflow regime and, a fortiori, in steady inflow regime, where a 
divergency occurs. For this reason, the effective yield cannot be considered 
alone, but together with the factor, ln[/x/(/x)j]. More specifically, the right- 
hand side of Eq. (12"51) . — p" ln[/x/(/x)j] = — pln[/i/(/i)j]/(l + k), is positive in 
all regimes due to the trend exhibited by the active gas mass fraction. In 
steady inflow regime, the following relation necessarily holds: 

_ Z Q - (Z )i _ t-tj 
{Zo)f — (Zo)i tf — ti 

in terms of the oxygen abundance related to a constant active gas mass 
fraction, or in terms of the cosmic time. 

In any case, mass conservation follows from integration of Eq. (129 p . as: 

fx + s + D = fi + s + D = fi = l ; (31) 
D = anS = ks ; (32) 

conformly to Eq. (|14p . where D is the gas mass fraction which outflowed from 
the box into the reservoir (k > 0) or inflowed into the box from the reservoir 
(k < 0). An equivalent form of Eq. (l3Tj) reads: 

li + s = 1 - D ■ (33) 

where variables on the left and the right-hand side member relate to the box 
and the reservoir, respectively. 

The birth-rate stellar function (number of stars born per unit volume, 
mass, and time), needs to be specified for determining the temporal behaviour 
of gas and star fractional masses. Following an earlier attempt (COO), the 
selected choice reads: 

B(m,t) = ^f(fz)$(m) ; m = — ; t= ; (34) 

V m Q Gyr 

where B is a normalization constant, V the volume of the box, $(m) the 
stellar initial mass function and f{p) an assigned function of the active gas 
mass fraction, which is the time-dependent term. The number of stars born 
per unit volume within an infinitesimal dimensionless mass range, m=F dm/2, 
and infinitesimal dimensionless time range, t =p dt/2, is B(m,t) dmdt. 

The number of long-lived stars generated (within the box) up to an as- 
signed dimensionless cosmic time, t, is: 

ft „ emu ~ _ 

N u (t)=Ldt l VB(m,t)dm = BF(t)I'(5) ; (35a) 

•Hi Jmjnf 

F(t)= Lf(fj)dt ; J'(5) = / $(m)dm ; (35b) 

Jti J m mi 



24 



where m m f is the lower stellar mass limit, and mu the upper mass limit of 
long-lived stars. 

The mass fraction in stars globally generated (within the box) up to a 
dimensionless cosmic time, t, is: 

S(t) = I dt I — B(m, t) dm = CF(t) ; (36a) 

Jt% J trimf Mq 

C = — -1(1) ; 7(1)= / m$(m)dm ; (36b) 

Mq Jmmi 

where m M f is the upper stellar mass limit. 

The mass fraction in long-lived stars and in gas outflowed from the box 
into the reservoir or inflowed into the box from the reservoir (in both cases 
with equal composition with respect to the preexisting gas), at an assigned 
dimensionless cosmic time, t, are: 

s(t) = aS(t) = aCF(t) ; (37) 
D(t) = anSit) = anCFit) ; (38) 

according to the definition of lock and flow, respectively. 

The star formation rate (within the box), the lock rate, and the flow rate 
at an assigned dimensionless cosmic time, t, are: 

= c/» ; ( 39 ) 

ds 

-^ = aCf(v) ; (40) 

^i = aK Cf(ri ; (41) 
dt 

and the combination of Eqs. ( 12 9 p and (I3"9~|) yields: 

a(l + K)CfQi) = -^ ; (42) 

which is equivalent to: 

r "AtL = - a (l + K )C(t-U) ; (43) 

in the special case, /(/i) = /i, Eq. (1431) can be integrated as: 

H = iaexp[-a(l + n)C{t-ti)} ; (44) 
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and the combination of Eqs. f[2T)j) . (|28p . and (|44p yields a linear trend for the 
oxygen abundance: 

Z -(Zo) i =paC(t-t i ) ; (45) 

regardless of the flow, k. 

In the more general case, /(/i) = ji v ', Eq. ( 13"§j) reduces to a Schmidt (1959, 
1963) star formation law. Leaving aside the above discussed special case, 
v — 1, Eq. (j4"3l can be integrated and, after some algebra, the result is: 



fJL = fJLi 



l-^z£a(l + K)C(t-U) 

Mi 



1/(1-") 



1/^1 



(46) 



which, using the logarithm Taylor series, ln(l + 2 /2 + x 3 /3- 

| x |< 1, and neglecting the terms of higher order with respect to the first, 
may be approximated as: 



H = ^ exp 



-Hi~ a(l + K)C(t - U 



1 - u)tf~ 1 a(l + n)C(t - U) <1 ; 



(47a) 
(47b) 



in the limit, v — > 1, Eqs. ( 144 p and (I47a0 coincide, as expected. The combina- 
tion of Eqs. (T20]), (F58D, and gZD yields: 



o)i 



PHiaCit - U 



(48) 



regardless of the flow, k. Then a Schmidt star formation law with exponent, 
v ~ 1, implies a linear dependence of the oxygen (or any other primary 
element) abundance on the cosmic time. 

In general, simple CBR models are described by Eqs. (I2"8"j) and (T431) . the 
latter only for including the temporal behaviour. 



3.2 Theoretical differential oxygen abundance distri- 
bution (TDOD) 

For a selected spectral class, the number of long-lived stars generated (within 
the box) up to an assigned dimensionless cosmic time, is: 

_ pt _ rm,2 

N(t) = I dt I VB(m,T) = C 12 F(t) ; (49a) 

Jti J mi 

C 12 = BI' (1,2) ; I'(l,2)=/ $(m)dm ; (49b) 

J mi 

where mi and mi are the lower and upper mass limit of the spectral class. 
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For simple CBR models, the oxygen abundance is monotonically increas- 
ing in time owing to Eqs. (I23p and ()46[) . which implies the number of stars 
born up to the cosmic time, t, coincides with the number of stars with 
oxygen abundance up to the related value, Zo(t), or N(t) = N(Zo). The 
same holds for the long-lived star mass fraction (including stellar remnants), 
s(t) = s(Z Q ). By use of Eq. (TT4"|) . the long-lived star mass fraction and the 
number of long-lived stars, N^, can be related as: 



mN u - m(N ei 



s f - si m(N u ) f - m(N u )i 



(50) 



where m is the mean mass of long-lived stars (including stellar remnants), 
fn = M s /Ni£. Under the assumption of a universal stellar initial mass func- 
tion, the fraction of long-lived stars belonging to a selected spectral class 
with respect to the total, maintains unchanged, N^/N = const, and Eq. (150]) 
translates into (Pagel and Patchett 1975): 



N-Nt 



sj-Si Nf-Ni 
and the combination of Eqs. fl3T|) and fl32|) yields: 

(1 + k)s — 1 — fx ; 
accordingly, Eq. f loT]) may be cast under the equivalent form: 

N-Nj _ 1 - n/fij 

N f -Ni~ l-fif/m ' 



(51) 



(52) 



(53) 



where the ratio, is expressed by Eq. (128|) with respect to the oxygen 

abundance, Z Q , and by Eqs. (jSJ) and (jlBjl with respect to the dimensionless 
cosmic time, t. 

In terms of the normalized oxygen abundance, = Zo/(Zo)o, Eq. f[28]) 
translates into: 



(54a) 
(54b) 



i = 




v" m 


E 




(Zoh 


Hi 


ii 


< 




< l 


Mi 









and the derivation of Eqs. (15T|) and (I5^|) with respect to yields: 
diV ds -d(/i//ii) 



(N f - Ni) 



(s f -Si)d(f) (1 - /j,f//ii)d(f> 



(55) 
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which, using Eq. (JMj) . takes the explicit form: 



dN = jij (Zq ) q 
(Nf — Ni) d(fi m-Hf p" 



exp 



(Zoh 



p" 



i) 



(56) 



where the decimal logarithm of the left-hand side is defined as the TDOD 
(Pagel 1989; COO; C01): 



ip((f)) = log 



dN 



(N f - N) d<f) 
1 (Zo'h _ 1 



acf) + b 
(Zo)e, 



log 



In 10 p" 

fit 



In 10 p 



1 + k) 



lnl0)a 



— at 



(57) 
(58) 
(59) 



which is represented by a straight line on the (O0?/>) plane. 

The TDOD slope, a, defined by Eq. ( |58|) . depends on the flow regime 
discussed above. More specifically, a < both in outflow regime and in 
moderate inflow regime; a = in steady inflow regime; a > in strong 
inflow regime. The TDOD intercept, b, defined by Eq. (l5"9~i) . must necessarily 
fulfill the condition, fif > 0, which implies the inequality: 



b > 6(/i/ = 0) = log(— In 10 a) — a(f>i ; a < 
on the other hand, fif > directly follows from a > 0. 



(60) 



3.3 Fitting to empirical differential oxygen abundance 
distribution (EDOD) 

Both the EDOD and the TDOD can be represented on the {0<pip) plane, 
but related normalizations are different. More specifically, the former is nor- 
malized to the sample population, N, according to Eq. (jHJ), while the latter 
is normalized to the computed long-lived star population, Nf — N, accord- 
ing to Eq. (j57p . Then the EDOD and the TDOD differ by a normalization 
constant, log Cn, which must be taken into consideration in performing the 
fitting procedure. In other words, the TDOD, ip^, has to be vertically shifted 
on the (00-0) plane by a value, logCV, for matching to the EDOD, ife. The 
TDOD intercept, defined by Eq. (1391 . translates into: 



b = lot 



lnl0)a 



a<pi ; 



(61) 
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where a and b can be determined as the EDOD regression line slope and 
intercept, respectively. 

Using Eq. (158ft . the empirical counterpart of Eq. (1271) . which implies the 
validity of Eq. ( 154"]) . reads: 



(^o)©l lnJ ^ 



In 10 a Hi 



< 



In 10 a fn 



< 1 



(62) 



keeping in mind that the active gas mass fraction, /i, monotonically changes 
as a function of the normalized oxygen abundance, 0, for simple CBR models. 

The particularization of Eq. ( I54p to the final configuration, by use of 
Eq. (|58|) . reads: 

L iln^ ; 



In 10 a Hi 

and the combination of Eqs. (16"T|) and (16"3"|) yields: 



i i 



In 10 a 



-exp lo (a0i + &){!- exp lo [a(0/ - (pi)]} ; 



(63) 



(64) 



which may be determined from the knowledge of the EDOD and related 
regression line. 



3.4 Different stages of evolution 

The mere existence of a G-dwarf problem in different regions of the Galaxy 
(e.g., van den Bergh 1962; Schmidt 1963; H76; Prantzos 2003; Ferreras et 
al. 2003) and perhaps in all galaxies (Worthey et al. 1996; Henry and 
Worthey 1999) implies the EDOD cannot be fitted by a straight line, as 
predicted by simple CBR models, but by a continuous broken line at most. 
To this aim, simple CBR models shall be extended by allowing different flow 
regimes during different stages of evolution, and defined as simple multistage 
closed- (box+reservoir) (MCBR) models. Accordingly, the flow, k, is different 
in different stages, while the equations of the model maintain their formal 
expression where variables and parameters are indexed by a letter, U = I, II, 
III, which denotes the stage under consideration. 

With regard to the U-th stage, the TDOD, defined by Eq. (|5T|) . reads: 

/ ix \ i (Cu)Ndsu 

K s u)/ - («u)ij d0u 

where (su)i an d ( s u)/ are the fractional star mass at the beginning and at 
the end, respectively, of the U-th stage, and log(C*u)N is the related nor- 
malization constant with respect to the EDOD, under the assumption of a 
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universal stellar initial mass function, according to Eq. (15 ip. The combina- 
tion of Eqs. d5TJ, §SD, and (JSSD yields: 



ifo{<fo) = «u0u + b v ; 

1 (Zq) q _ 1 (Z ) e 



au = 

b v = log 



In 10 p'{j In 10 p v 

(C<u)N(/iu)l (-lnlO)au 



_(/iu)i - (/X U ) / 
which is valid provided the following inequality holds: 



In 10 au (/iu) 

in agreement with Eq. (162 p . 
The particularization of Eq. 

(<hh - 



< 



( z o)e i ln Mf 



In 10 au (/Ux 



< 1 



to the U-th stage reads: 



1 1 ln (Ri)/ 



In 10 au (/iu)i 
and the combination of Eqs. (167]) . f )68|) . and fTTOj) yields: 
1 1 



(66) 
(67) 

(68) 



(69) 



(70) 



(G 



UjN 



In 10 au 



exp lo [a u (0u)i+M{ 1 - ex Pio{ a u[(0u)/-(0u)i]}} ; (71) 



which may be determined from the knowledge of the EDOD belonging to the 
U-th stage and related regression line. In general, different samples and/or 
different stages imply different values of the normalization constant. 



3.5 Application to a special stellar system 

For selected [0/H]-[Fe/H] relations, the EDOD related to the fslO sample 
(taken as representative of ther inner Galactic halo) can be divided into four 
or three regions where the trend is linear to a good extent, as shown in 
Figs.H] and [5j In the light of the model, region A corresponds to the first 
stage of evolution, where the system is still assembling, characterized by 
strong inflow regime (k < —1, a > 0); region F corresponds to the second 
stage of evolution, where the system is forming, characterized by outflow 
regime (k > 0, a < 0); region C corresponds to the third stage of evolution, 
where the system is undergoing contraction, characterized by outflow regime; 
region E corresponds to the fourth step of evolution, where the system has 
attained dynamical equilibrium, characterized by outflow regime. 
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The gas and star mass fraction are left unchanged passing from a selected 
stage to the next one, which implies the validity of the following relations: 



(Xv)f = (Xv)i ; (X A )i = Xi ; (X E ) f = X f ; (72a) 
X = fi,s,D ; U = A,F,C ; V = F, C, E ; (72b) 

where X iy Xf, are related to the whole evolution regardless of the stages. 
Accordingly, mass conservation during the U-th stage may be expressed as: 

/iu + (1 + k v )s v = fly + s v + Aj = (/iu)i + (su)i + (-Du)i ; (73) 

where the initial values are known via Eq. (1721) . 

The regression line in an assigned region of the EDOD is defined by the 
slope, ay, and the intercept, b\j, and the intersections between regression 
lines related to adjacent regions mark initial and final values of normalized 
oxygen abundance, {4>v)i an d (4>u)f- The normalization constant, (Cu)n> 
is determined via Eq. (T7T]) . The final value of the active gas mass fraction, 
(/iu)/, follows from Eq. (T70]) . 

The yield, p, and the lock, a, depend on the birth-rate stellar function. 
For a power-law stellar initial mass function, the EDOD in different Galactic 
environments is reproduced to a good extent from simple CBR models with 
the following output parameters (CM09): 

1.0340 ; (74a) 



(Zoh 

a 2 . 9 = 0.66604 ; a 2 . 35 = 0.85360 ; (74b) 
(m mf ) 2 . 9 = 0.44449 ; (m mf ) 2 . 35 = 0.015436 ; (74c) 

for an assumed solar oxygen abundance, (Z ) Q = 0.0056, where the indices 
relate to power-law exponents equal to 2.9 and 2.35, respectively. For further 
details refer to earlier attempts (C01; C07). 

The flow, k\j, is deduced from Eq. (1671) . The final star mass fraction, 
(s\j)f, and outflowed or inflowed gas mass fraction, (-Du)/, are determined 
from Eq. (1731) particularized to the end of the U-th stage. 

In the light of the model, the initial star mass fraction, s», results from the 
fractional mass in stars (including stellar remnants) with normalized oxygen 
abundance, <p < (pi. The above value, though still uncertain at present, may 
be understimated to a good extent as Sj = 0. In addition, Di = without 
loss of generality, which implies /Xj = 1. 

The fractional mass of the box (with respect to the initial value) attains 
the maximum value at the end of A stage, where gas inflows into the box 
from the reservoir. The related value is (//a)/ + ( s a)/- The fractional mass 
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in stars at the end of evolution, which coincides with the end of E stage, is 
(se)/- Accordingly, the mass ratio of the box at the end of evolution to the 
outflowed gas reads: 



provided the earlier stage A is in inflow regime and the subsequent stages 
F, C, (or FC), E, are in outflow regime, as it is the case for the EDOD and 
related TDOD under consideration. 

With the above values of input parameters, the flow, the active gas mass 
fraction, the star mass fraction, and the outflowed or inflowed gas mass frac- 
tion, at the end of each stage of evolution, can be computed. The results are 
listed in Table for the EDOD related to the fslO sample for [0/H]-[Fe/H] 
empirical relations deduced from the RB09 sample (top panel) and Fal09 
sample, case SHI (middle and bottom panel), see Fig.|H With regard to the 
bottom panel, stages F and C are merged into a single one, FC. 

In the former case (top panel), where the empirical [0/H]-[Fe/H] relation 
has been determined using the LTE approximation (RB09), the mass of the 
box is increased by a factor of about 10 at the end of A stage, and is reduced 
to about the initial value at the end of E stage. The mass ratio, Mb ox /M fl, 
is about one tenth via Eq. (175|) . 

In the latter case (middle and bottom panel), where the empirical [O/H]- 
[Fe/H] relation has been determined relaxing the LTE approximation (Fal09), 
the mass of the box is increased by a factor of about 3-4 at the end of A 
stage, and is reduced to about one half the initial value at the end of E stage. 
The mass ratio, M^ ox /M fi, is about one sixth via Eq. (1T51) . 

In any case, the initial mass of the box at the beginning of A stage (proto- 
inner halo) is comparable to the mass of the box at the end of E stage (present 
inner halo). 

3.6 Discussion 

Though TDODs calculated using simple MCBR models of chemical evolution 
fit to EDODs with different linear trends in different regions, as shown in 
Fig.EJ still the application to the inner Galactic halo remains speculative in 
absence of further improvement. The main reasons are outlined below. 

First, the EDOD has been derived from a fictitious sample, fslO, which, 
in turn, has been deduced from two incomplete samples, RN91 and H V, both 
biased (to a different extent) towards low metallicities and affected by disk 
contamination. A homogeneous and unbiased sample would be needed to 



M box 
M ofl 



M/ + (s E )/ 



(75) 



(>a)/ + (s A )f - M/ - (s E )/ 
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Table 7: Input parameters (deduced from the regression lines), (0u)i> a u> &u, 
and output parameters, k\j, (nv)f, ( s u)/> (Dv)f, for simple MCBR models 
fitting to the EDOD related to the fslO sample, for [0/H]-[Fe/H] empirical 
relations deduced from the RB09 sample (top panel) and Fal09 sample, case 
SHI (middle and bottom panel). Four (A, F, C, E; top and middle panels) or 
three (A, FC, E; bottom panel) stages of evolution are considered, according 
to the linear trends exhibited by the EDOD (Fig.HJ). Stages O before A and 
after E are not considered as no sample object lies within the corresponding 
metallicity range. Other input parameters are (Zq)q = 0.0056; p/(Zo)® = 
1.0340; Hi = (/io)/ = 1; Si = (so)/ = 0; Di = (-Do)/ = 0; where the index, O, 
denotes oxygen with regard to the solar abundance and stage O of evolution 
otherwise. For further details refer to the text. 



U 



5U i 



flu 



K V 



(/i 



( g u)/ 



A 


9.4624E- 


-4 


+1.1382E+2 


-1.0567E- 


1 


-2.7201E+2 


9.9175E- 


-0 


3.2905E- 


-2 


F 


9.7001E- 


-3 


-2.0950E+0 


+1.0188E- 





+3.9881E+0 


7.7178E- 


-0 


4.7391E- 


-1 


C 


6.1687E- 


-2 


-7.3565E+0 


+1.3433E- 





+1.6515E+1 


4.5710E- 


-1 


8.8844E- 


-1 


E 


2.2854E- 


-1 


-2.2569E+0 


+ 1.7784E- 


1 


+4.3735E+0 


3.5961E- 


-2 


9.6681E- 


-1 


O 


7.1779E- 


-1 


















A 


3.1623E- 


4 


+7.9432E+1 


+2.7571E- 


1 


-1.9012E+2 


3.6714E- 


-0 


1.4125E- 


■2 


F 


7.4271E- 


3 


-3.8855E+0 


+8.9453E- 


1 


+8.2913E+0 


1.2551E- 


-0 


2.7531E- 


1 


C 


1.2740E- 


1 


-2.6460E+0 


+7.3661E- 


1 


+5.3000E+0 


5.4973E- 


-2 


4.6581E- 


-1 


E 


6.4083E- 


1 


-5.3661E-1 


-6.1514E- 


1 


+2.7764E-1 


2.4384E- 


-3 


5.0693E- 


1 


O 


3.1623E- 





















A 


3.1623E- 


-4 


+7.9432E+1 


+2.5751E- 


1 


-1.9012E+2 


3.2209E- 


-0 


1.1743E- 


-2 


FC 


6.7114E- 


-3 


-2.8643E+0 


+8.2804E- 


1 


+5.8196E+0 


5.6405E- 


-2 


4.7577E- 


-1 


E 


6.2002E- 


-1 


-5.3661E-1 


-6.1514E- 


1 


+4.7764E-1 


2.4384E- 


-3 


5.1801E- 


-1 


O 


3.1623E- 


-0 
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this respect. The EDOD high-metallicity tail could be largely due to disk 
contamination as stage E begins at [Fe/H] ~ —1, which is close to both the 
transition from halo to bulge/disk globular cluster (e.g., Mackey and van den 
Bergh 2005) and the peak of bulge EDOD (Sadler et al. 1996). Different 
linear trends in different regions are also shown by globular clusters (C07), 
bulge (C07), and disk (Caimmi 2008) EDOD. 

A linear trend exhibited by the EDOD related to different Galactic en- 
vironments or galactic environments provided the G-dwarf problem is uni- 
versal (Worthey et al. 1996; Henry and Worthey 1999), implies outflowing 
or inflowing gas with same composition as the preexisting gas. If the above 
mentioned effect is real, the chemical evolution within the building blocks 
of galactic subsystems, in particular the inner Galactic halo, occurred to a 
similar extent. 

Second, using different methods yields different [0/H]-[Fe/H] empirical 
relations, as shown in Fig.[3l even in dealing with coinciding sample objects, 
as listed in Table [TJ With regard to the inner Galactic halo, the initial 
mass deduced using two selected [0/H]-[Fe/H] empirical relations, has been 
found to change by a factor of about two. Future improvement would imply 
direct oxygen abundance determinations (Ramirez et al. 2007) where differ- 
ent methods yield consistent results or, at least, a [0/H]-[Fe/H] empirical 
relation which attains general consensus. 

Third, though the assumption of null star mass fraction holds to a good 
extent for the starting configuration, still the presence of stars with lower 
metallicity with respect to sample objects, [Fe/H] < —4.2, has to be consid- 
ered. This extremely low metallicity has currently been detected in about 
a dozen of stars down to [Fe/H] « —5.4 (e.g., Beers and Christlieb 2005). 
Basing on theoretical arguments, a lower limit oxygen abundance in Pop II 
stars has been determined as [O/H] = — 3.05 =F 0.20 (Bromm and Loeb 2003), 
where earlier nucleosynthesis comes from more massive Pop III stars. 

If the transition from the latter to the former population was not instan- 
taneous, and coeval Pop III and Pop II stars were generated (Smith et al. 
2009), then the stellar initial mass function can no longer be considered as 
universal in time for Pop II stars. Accordingly, model evolution must be 
started after the last Pop III star has undergone supernova explosion, which 
has been assumed to take place when the metal abundance is [Fe/H] = —4.2. 
In absence (to the knowledge of the author) of a reliable estimate of the 
mass fraction in stars and stellar remnants with initial metal abundance, 
[Fe/H] < —4.2, a null value has been assigned. 

For a stellar system resembling the inner Galactic halo, a nontrivial ques- 
tion is concerned with the nature of the reservoir where gas outflows after 
the assembling (A) stage. A natural candidate is the bulge, as suggested by 
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comparison between specific angular momentum distributions related to halo 
and bulge stars (Wyse and Gilmore 1992). According to recent attempts, (1) 
the Galactic halo is made of an inner component which exhibits a modest 
prograde rotation and a metallicity peak at [Fe/H] ~ —1.6, and an outer 
component which exhibits a clear retrograde net rotation and a metallicity 
peak at [Fe/H] rj -2.2 (Carollo et al. 2007, 2010); and (2) the Galactic 
bulge is made of an inner component which exhibits bar-like kinematics and 
metal-rich population, and an outer component which exhibits spheroid kine- 
matics and metal-poor population (Babusiaux et al. 2010). Then a similar 
specific angular momentum distribution related to inner halo and outer bulge 
is expected. 

In this view, the mass ratio of the box at the end of evolution to the gas 
outflowed from the box into the reservoir during the evolution, Mb ox /M fl, 
expressed by Eq. (175]) . should reproduce the mass ratio of (inner) halo to 
(outer) bulge, as: 

(M H )inn = (MhWMh M h 

(M B )out (MbWMbMb ' 1 J 

which is equivalent to: 

(M H )inn = (M R ) inn /(M B ) out (Mgjout . . 

M H Mr/Mb M b [ ' 

where the indices, inn and out, denote the inner and the outer compo- 
nent, respectively, H the halo, and B the bulge. The the inner halo to 
the outer bulge mass ratio, deduced from the results listed in Table [7J is 
(M H )inn/(MB)out ~ 0.10-0.20. The halo to bulge mass ratio is currently esti- 
mated as Mh/Mb ~ 0.05-0.10. Then simple MCBR models considered here 
yield an inner halo fractional mass (normalized to the halo) which is compa- 
rable with, or exceeding by a factor up to 4, the outer bulge fractional mass 
(normalized to the bulge). On the other hand, quantitative results cannot 
be expected for the above mentioned reasons. In conclusion, the (inner) halo 
to (outer) bulge mass ratio appears to be an additional output parameter 
provided by simple MCBR models of chemical evolution, with regard to the 
inner Galactic halo. 

In the light of the model, the inner Galactic halo can be considered in- 
dependently of the outer halo and the bulge for the following reasons. The 
outer halo is both dynamically and chemically decoupled from the inner halo 
(Carollo et al. 2010), which implies accretion by smaller subunits (mainly 
dwarf spheroidal galaxies) after the formation of the inner halo. In fact, 
globular clusters belonging to the inner and outer halo are usually classi- 
fied as "old halo" and "young halo", respectively (e.g., Mackey and van den 
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Bergh 2005). Similarly, a (inner) halo - (metal-poor) bulge collapse can be 
inferred by comparison of related specific angular momentum distributions 
(Wise and Gilmore 1992), which implies bulge formation is subsequent to 
inner halo formation, and then the (metal-poor) bulge may be considered as 
a reservoir for the outflowing gas, in the sense specified by the model. 

4 Conclusion 

Under the assumption that two samples of halo stars, RN91 and HV, are 
equally representative of the inner Galactic halo within the metallicity range, 
—3.0 < [Fe/H] < —2.8, a fictitious sample, fslO, has been built up and taken 
as representative of the inner Galactic halo within the metallicity range, 
—4.2 < [Fe/H] < +0.2. The related differential empirical oxygen abundance 
distribution (EDOD) has been established using different [0/H]-[Fe/H] em- 
pirical relations, deduced from different samples (RB09, Fal09, Sal09) where 
different methods have been exploited for determining the oxygen abundance. 

More precisely, the EDOD has been deduced from the fslO sample by use 
of two alternative [0/H]-[Fe/H] empirical relations: one, determined from 
the RB09 sample in presence of local thermodynamical equilibrium (LTE) 
approximation (RB09), and one other, determined from the Fal09 sample 
in absence of LTE approximation with due account taken of the inelastic 
collisions via neutral H atoms, SHI (Fal09, Sfi = 1), as shown in Fig. [31 

A linear trend has been exhibited by related EDODs within three or 
four regions, as shown in Fig.Hl and the slope and intercept estimators of 
corresponding regression lines have been determined together with their dis- 
persions, as listed in Table [2] using different interpolation methods. It has 
been pointed out that the earlier trend, characterized by positive slope, is a 
signature of a G-dwarf problem. 

The main uncertainties on the EDOD have been recognized as related 
to (1) biases on the RN91 and HV samples due to selection effects towards 
sufficiently low ([Fe/H] < -4.2) and sufficiently high ([Fe/H] > -2.8) metal 
abundance (HV) and disk contamination for [Fe/H] > —2.0 (RN91; HV), 
and (2) lack of clear indications on a recommended method for determining 
oxygen abundance, as shown in Table [Hand in Fig. [1] for 11 stars in common 
among the RB09 and the Fal09 sample. 

Fitting a theoretical differential oxygen abundance distribution (TDOD) 
to the above discussed EDOD has needed an extension of simple closed-box 
(CB) models of chemical evolution on two respects. First, the system has 
been conceived as made of a box and a reservoir, where the following pro- 
cesses have been allowed: gas outflow from the box into the reservoir (H76), 
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moderate gas inflow (C07), steady and strong gas inflow (current paper) into 
the box from the reservoir. Simple closed- (box+reservoir) (CBR) models 
have exhibited mass conservation within the system (box+reservoir), while 
it has been violated within a single subsystem (box or reservoir). Second, the 
history of the system has been conceived as a succession of different stages 
characterized by different outflow or inflow rate. Simple multistage closed- 
(box+reservoir) (MCBR) models have been found to yield TDODs in the 
form of continuous broken lines, which can fit to related EDODs. 

An application of MCBR models to a special stellar system resembling 
the inner Galactic halo has been made with fiducial values of input param- 
eters which cannot be deduced from the EDOD. The metal abundance at 
the beginning and at the end of each stage, have been inferred from the 
intersection of regression lines fitting to adjacent regions of the EDOD, as 
shown in Fig.[5j The mass ratio of the box at the end of evolution to the gas 
outflowed from the box into the reservoir through the evolution, has been 
determined as M box /M ofl m 0.11 for a [0/H]-[Fe/H] empirical relation de- 
duced from the RB09 sample in presence of LTE approximation (RB09), and 
Mbox/MoQ ~ 0.16-0.19 for a [0/H]-[Fe/H] empirical relation deduced from 
the Fal09 sample in absence of LTE approximation with due account taken 
of the inelastic collisions via neutral H atoms, SHI (Fal09, Sfi = 1), with 
regard to four and three stages of evolution, respectively. 

For current estimates of the halo-to-bulge mass ratio, M ha i /M bu i ge « 
0.05-0.10, the inner halo fractional mass (normalized to the halo) has been 
found to be comparable with, or exceeding by a factor up to 4, the metal- 
poor bulge fractional mass (normalized to the bulge). On the other hand, 
it has been considered that quantitative predictions cannot be made for the 
Galaxy unless complete and unbiased samples of the inner Galactic halo are 
available, and discrepancies among [0/H]-[Fe/H] empirical relations related 
to different samples and different methods are removed. 
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Appendix 
A Tables 

Tables [31 HI and [7] may be exceedingly large in the text, and for this reason 
they are broken to be completely accessible. 
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Table 8: Table El left. 



[Q/H] = 0.72 [Fe/H] (RB09) 
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Table 9: Table El right. 
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Table 10: Table H left. 
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Table 11: Table H right. 
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Table 12: Tabled left. 
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Table 13: Tabled right. 
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Figure 1: Comparison between (i) [O/H]- [Fe/H] empirical relations deduced 
from 11 halo stars in common among RB09 (upper panels, diamonds) and 
Fal09, case LTE (top left, squares), case SHI (top right, crosses), and (ii) 
[Fe/H] (bottom left) and [O/H] (bottom right), case LTE (squares), case 
SHI (crosses), deduced from the above mentioned stars. The composite 
symbols mark LP831-070, where only an upper limit to oxygen abundance is 
available in the Fal09 sample. The straight lines in upper panels are [O/H] = 
0.72[Fe/H] (left) and [O/H] = [Fe/H] +0.70 (right). The 1:1 correspondence 
in lower panels is represented by a dashed line. Typical uncertainties are 
visualized as crosses on the top left corner of each panel. 
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Figure 2: [Fe/H] distributions (normalized to the unit area) related to four 
different samples discussed in the text. Sample denomination and population 
are indicated on each panel. In all cases, the bin width is A [Fe/H] = 0.2. 
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Figure 3: The [0/H]-[Fe/H] empirical relation deduced from two samples 
with heteroscedastic data, RB09 and Sal09, and three samples with ho- 
moscedastic data, Fal09, cases LTE, SHO, and SHI, indicated on each panel 
together with related population. Also plotted on each panel are the adopted 
linear dependences, [O/H] = 0.72 [Fe/H] (full) and [O/H] = [Fe/H] + 0.70 
(dashed). The regression lines related to five different methods are shown for 
each sample on the top right panel. For further details refer to the text. 
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Figure 4: The empirical, differential oxygen abundance distribution (EDOD) 
related to the fslO sample for [0/H]-[Fe/H] empirical relations deduced from 
the RB09 (left panels) and Fal09, case SHI (right panels) sample. The whole 
distribution is represented in upper panels, while lower panels zoom the low- 
metallicity range. The uncertainty of the distribution is determined from 
Poisson errors. The vertical bars on the horisontal axis mark [Fe/H] = —2.2, 
— 1.6, and —0.8. The vertical dotted lines mark [Fe/H] = —2.7, —1.7, and 
—0.9, where the linear trend of the EDOD changes passing from a region to 
the adjacent one. For further details refer to the text. 
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Figure 5: Regression lines to the empirical differential oxygen abundance 
distribution (EDOD) plotted in Fig.HJ with regard to the regions (from the 
left to the right): A (dotted, positive slope, tia = 6), F(dotted, negative 
slope, n F = 6), C (full, n c = 5), E (dashed, tie = 4). The more inclined 
dashed line on right panels corresponds to the FC (npc = 10) region, in 
alternative to F and C regions separately. Other captions as in Fig.|H For 
further details refer to the text. 
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